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ABSTRACT 

We report the discovery of a UV-bright tidal dwarf galaxy candidate in the NGC 4631/4656 galaxy 
group, which we designate NGC 4656UV. Using survey and archival data spanning from 1.4 GHz to 
the ultraviolet we investigate the gas kinematics and stellar properties of this system. The H I mor- 
phologies of NGC 4656UV and its parent galaxy NGC 4656 are extremely disturbed, with significant 
amounts of counterrotating and extraplanar gas. From UV-FIR photometry, computed using a new 
method to correct for surface gradients on faint objects, we find that NGC 4656UV has no significant 
dust opacity and a blue spectral energy distribution. We compute a star formation rate of 0.027 Mq 
yr~^ from the FUV fiux and measure a total H I mass of 3.8x10* Mq for the object. Evolutionary 
synthesis modeling indicates that NGC 4656UV is a low metallicity system whose only major burst of 
star formation occurred within the last ^ 260 — 290 Myr. The age of the stellar population is consis- 
tent with a rough timescale for a recent tidal interaction between NGC 4656 and NGC 4631, although 
we discuss the true nature of the object-whether it is tidal or pre-existing in origin-in the context 
of its metallicity being a factor of ten lower than its parent galaxy. We estimate that NGC 4656UV 
is either marginally bound or unbound. If bound, it contains relatively low amounts of dark matter. 
The abundance of archival data allows for a deeper investigation into this dynamic system than is 
currently possible for most TDG candidates. 

Subject headings: galaxies: individual (NGC 4656) - galaxies: dwarf- galaxies: interactions - galaxies: 
kinematics and dynamics - galaxies: structure 



1. INTRODUCTION 

Systems of interacting galaxies have long been known 
to contain some very irregular structur es such as rings, 
bridges, plumes, and tidal tails (see lSchweize"il 119861 
for a review). In recent years, tidal dwarf galaxies 
(TDGs) have also been recognized as byproducts of col- 
lisions between massive galaxies, identified by regions 
of active s tar formation a nd res ervoirs of cold gas (see 
iDuc fc Mirabel 1999: d1I3 120111 for a review). Simula- 
tions show that they consist of baryonic material stripped 
fro m the outer disks of interacting gas-rich 'p arent' galax- 
ies (IHibbard fc Mihosl[T995t fPuc et al.l[2000l) . TDGs are 
believed to form stars through gas instabilities and col- 
lapsi ng H I clouds, and contain little to no dark mat - 
ter (jBarnes fc HernquistI Il992t lElmegreen et all 119931 ). 
In general, TDGs appear to fall into one of two cat- 
egories based on their stellar population: those domi- 
nated by old stars ripped from their parent galaxies, and 
those that consist almost exclusively of new stars forming 
from the H I gas that wa s stripped during the in terac- 
tion (Due fc Mirabe]|ll9'99l : iDuc et al.i 2000 : Braine et al.l 
[2OOII) . 

Tidal features are interesting for several reasons. They 
revea l the history of galaxy interactions (IHibbard et al.l 
Il994f ) and they may be used to trace the underly- 
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ing d ark matter distrib ution of the group or cluster 
(e.g. iMihoset al.]|1998[ ). The ultimate fate of TDGs 
is unknown, but due to their dearth of dark matter 
they should be minimally bound by self-gravity and 
therefore easily disrupted by further galaxy encoun- 
ters. They provide a potential source for metal en- 
richm ent of the intergalactic medium (jWiersma et al.l 
l2010f ) as the material they contribute will have a sim- 
ilar metallicity t o the i r alre ady-evolved parent galaxy 
puc fc Mirabell 11994 Il999f l. T hey are also a pro- 
posed origin of in tracluster light (jGregg fc WestI 119981 : 
iMihos et all 120051) . Finally, TDGs provide a laboratory 
for studies of star formation in low density environments 
which ar e not complicated b y the dynamics of spiral 
galaxies (|Boquien et al.l[2009| ). 

Neutral hydrogen observations play a crucial role in 
identifying TDGs. The H I gas in late-type galax- 
ies typically extends far beyo nd the stellar disk (e.g. 
Ivan der Kruit fc Shostakiri984f) . so it is easily disturbed 
or removed in gravitational encounters, revealing the 
kinematics of an interaction. Numerical simulations 
have successfully reproduced the distribution of gaseous 
and stellar tidal debris in interacting systems (e.g. 
iToomre fc Toomr^ll972D . A larg e fraction of this ga s 
may fall back to the parent galaxy (jHibbard et al.|[T994) , 
but some portion of it is pulled out with stars. As a re- 
sult, we observe gaseous debris among tidal streams, and 
TDGs typically contain their own reservoir of H I gas 
out of which they are currently forming stars. Instabil- 
ities in this gaseous compone nt may be the driving fac 
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tor in the formation of TDGs (iBarnes fc Hernguis.. __ . 
lElmegreen et al.ll993l: IDuc fc Brinksl200lD . The velocity 
dispersion of H I may be used to estimate the dynamical 
mass and evaluate whether H I density enhancements are 



gravitationally bound. However, the H I can be entan- 
gled in tidal material. CO gas, which forms in-situ, may 
provide a m ore accurate estimate of the TDG structure 
IjBraine et al. 2001). 

Generally, TDGs are observed to be young systems 
which may be an ind icator of their short lifetime 
(jHibbard fc MihosI IT995[) or ma y be due to th e diffi- 
culty in identifying older TDGs |Duc et al.ll2"007l ). How- 
ever, they are also generally observed to be chemically 
evolved: TDGs contain higher metallicities than normal 
dwarf galaxies, up to approximately | solar, and they are 
abund ant in molecular ga s as measured by CO observa- 
tions (jBraine et al.|[200ll ). Thus, metallicity estimates 
are used to confirm that TDG material is similar to that 
of its putative parent galaxy while population synthesis 
modeling is used to show that the stellar ages are less 
than the age of the ti dal feature to which they belong 
(e.g. iWeilbacher eTall 12000). 

How star formation proceeds in TDGs remains un- 
known and may provide valuable clues f o r how star for- 
mation occurs in general. iBoquien et al.l ()2009() find evi- 
dence that the star formation properties in intergalactic 
star forming regions are similar to that in spiral galax- 
ies. However, other observations suggest that star for- 
mation in low density environments, such as in dwarf 
galaxies or the extreme outer disks of massive galaxies, 
may take place below classical star formation gas sur- 
face density thresholds put forth by Kcnnicutt (1998) 
and Martin & Kcnnicutt (2001) (e.g. IThilker et al.1120071: 
iRovchowdhurv et al. 2009 ). Further, scatter in the pre- 
dicted star formation rates from UV, Ha, and 8.0 /xm 
observations in these low density environments is not un- 
derstood, but may be a function of age, metallicity, gas 
content, dust content, or the stellar initial mass function 
(JThilkcr 2011). As dynamically simple systems with a 
range of metallicities, molecular gas, and dust content, 
dwarf and tidal dwarf galaxies may be good laboratories 
for shedding light on the mechanics of star formation. 

We serendipitously discovered a UV-bright intergalac- 
tic star forming region, which was previously un- 
recognized or unacknowledged in optical observations, 
while examining GALEX images of the edge-on galaxy 
NGC 4656. We consider the object (hereafter referred 
to as NGC 4656UV) a strong candidate for tidal dwarf 
galaxy status. 

NGC 4656UV is a member of the loose (total number 
of members ~6; [Garcia 1993) galaxy group consisting 
primarily of two edge-on spiral galaxies, NGC 4656 and 
NGC 4631, which are suspected to have had a strong 
interaction in the past. An H I bridge discovered by 
[Roberts! (|1968f ) extends between the two spirals suggest- 
ing gas has been stripped from NGC 4631 by its gravita- 
tional encounter with NGC 4656. High resolution inter- 
ferometric observations of the group resolve the bridge 
into multiple H I "fi ngers" (iWehachew et al.| [T978: Rand 
[1991 . In particular, [Randl ([199D notes that NGC 4656 
contains evidence of a disturbed morphology resulting 
from a prior interaction, including non-circular motions 
and asymmetry. In the optical, previous attention has 
been paid to the sharp bend at the northeast edge of 
NGC 4656, a feature which gives the galax y its d istinct 
"hockey stick" appearance. [Stavton et al.[ ([1983[ ) used 
photometry to determine that the bend is a feature of 
the galaxy resulting from an interaction, as opposed to 



a background galaxy (NGC 4657). 

An ultra- faint Ha "sheet" has been detected b etween 
NGC 4656 and NGC 4631 ([Donahue et al.[[l995l) . how- 
ever NGC 4656UV is not detected in these observa tions. 
H I total intensity maps presented by I Randl ([1994*) show 
a density enhancement at the location of NGC 4656UV, 
however they suggest it is part of an H I ring surrounding 
NGC 4656 and do not identify the gas as being distinct 
from the parent galaxy. 

While nearly invisible in the Digitized Sky Survey 
and very close to detection limit s in the r ed {i,z) 
Sloan Digital Sky Survey (SDSS, lYork et al.l ([1)00)) 
bands, GALEX observations easily reveal the TDG can- 
didate. NGC 4656UV is located just to the northeast 
of NGC 4656 at 12^44"14.794'' -F32°16™48.26" (J2000) 
with a center-to-center separation of ^ 11 kpc, assuming 
a distance of 7.2 Mpc ([Seth_et al.'^2005). NGC 4656 and 
NGC 4656UV are connected by a faint bridge of H I gas 
and stars. At this distance NGC 4656UV is one of the 
most nearby examples of vigorous star formation occur- 
ring in an extreme, low density environment. 

In this paper we present our investigation into the 
nature of this UV-bright tidal dwarf galaxy candidate 
through the accumulation of publicly available survey 
and archival data. We make the case that it is most likely 
a tidal dwarf galaxy created from an encounter between 
NGC 4656 and NGC 4631, instead of a pre-existing, but 
highly rejuvenated, dwarf galaxy, or a tidal stream. The 
paper is organized as follows: we describe the archival 
data and our post-processing methods in ij2l we present 
our results in U3|and discuss their implications in 311 We 
end with concluding remarks in ij5| 

2. DATA 

A plethora of survey and archival data exists on 
NGC 4656 and its companion candidate TDG across 
a broad wavelength regime, giving us nearly complete 
wavelength coverage from 2.01 keV to 21 cm. In 
addition to the discovery GALEX tiles, we obtained 
images f rom the SDSS, Two Micron All-Sky Survey 
(2MASS, ISkriitskie et all EUM )). Infrared Array Cam- 
era (IRAC) and Multi-band Imaging Photometer (MIPS) 
on the Spitzer Space Telescope, and archival H I obser- 
vations from the Very Large Array (VLA). ROSAT X- 
ray imaging of NGC 4656 includes the TDG candidate, 
however there is no hard or soft X-ray emission from 
NGC 4656UV down to the detection limi ts of the ob- 
serva tion ('^ 9 X 10^^ erg s~^; Figure 8 of iVogler et al.l 
[1991) . 

2.1. HI 21 cm observations 

We retrieved archival VLA data from the online 
National Radio Astronomy Observatory Science Data 
Archive^. The observations consist of a single pointing 
in D-array taken on 18 December 1989 (Project Code 
AR0215). It is centered on the optically bright portion 
of the edge-on spiral galaxy, NGC 4656 (12'M3'"57.536'' 
-f 32°10'06.5", J2000). The primary beam of the VLA at 
L-band is half a degree in diameter, such that NGC 4656 
and NGC 4656UV span nearly the full extent of the 
beam. The galaxy was observed for approximately six 

■^ https://archivc.nrao.edu/archivc/advquery.jsp 
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Fig. 1.— GALEX FUV image of NGC 4656 and NGC 4656UV 
with the apertures used for photometry indicated. The smaller el- 
lipse is for NGC 4656UV while the larger ellipse denotes NGC 4656. 
The boxes show the regions used for (from North-East to South- 
West) the Hook, Main, and Tail areas of NGC 4656. 

hours on source with a velocity resolution of 10.3 km s~^ . 
We calibrated and imaged the data in AIPS using the 
standard procedures. The resulting data cube has an 
RMS noise of 0.81 mJy beam~^ channel"-'^. At 5cr, we 
achieved a mass sensitivity of 5.10 x 10^ M© channel"^ 
and a column density of 2.51 x 10^^ cm~^ channel^^. 

The NGC 4631/NGC 4656 group of ga l axies has been 
previously o bserved by iWeliachew et al.l (|1978D and by 
iRandl ()1994[ ) , both with multiple pointings of the West- 
erbork Synthesis Radio Telescope (WSRT). The archival 
VLA observations presented here are a single pointing fo- 
cused on NGC 4656 and offer similar spectral and spatial 
resolution to the WSRT observations. However, the VLA 
observations provide an improve ment of 1.7 and 1. 4 times 
the sensitivity co mpared to the IWeliachew et al.l (|1978[ ) 
and IRandl ()1994| ) WSRT observations, respectively. 

2.2. UV, optical, & mid-IR imaging 

We visually inspected all UV, optical, and infrared im- 
ages to ensure that they were free from defects and con- 
tained both the main galaxy and the TDG candidate 
object. Only about half of NGC 4656UV is on Spitzer 
IRAC 4.5 and 8.0 ^m images, so we exclude those bands 
from photometry of the TDG candidate. To ensure com- 
plete and uniform coverage of the objects we downloaded 
image mosaics of the SDSS and 2MASS data using the 
Montage mosaicking software. 

An initial list of sourc es were extracted usin g the 
Source-Extractor package (jBertin fc Arnoutill996D . Be- 
cause SExtractor is not optimized for nearby extended 
sources, we hand-edited catalogs to remove spuri- 
ous/overlapping ellipses and add masking regions over 
unblanked stars and background galaxies. To ensure uni- 
form processing of the images, the hand-edited catalogs 
were homogenized for GALEX, SDSS, and MIPS bands. 
Due to the high number of contaminant point sources in 
the IRAC frames, additional masking was done for those 
bands. These extra sources were not blanked in the other 
images to maximize the accuracy of the photometry. In 
all cases we chose to mask sources where we were unsure 
if they were contaminants. 

We photometered both NGC 4656 and NGC 4656UV. 
We set elliptical apertures around both objects, with 
sizes chosen to encompass all UV emission. We also 



photometered three separate regions of NGC 4656, rep- 
resenting visually distinct areas of the galaxy. On the 
northeast side of NGC 4656 (toward NGC 4656UV) the 
'hook' is a significant 90° warp seen in images from the 
UV to the far-infrared. This feature, part of which is also 
known as NGC 4657, appears very dense and blue. We 
separate the central portion of NGC 4656 as the 'main' 
body. While there is no obvious bulge in NGC 4656 to 
denote the center of the system, this region has both the 
reddest optical colors and the highest surface-brightness 
of the galaxy. Finally, we photometer the extended 'tail' 
on the Southwest side of NGC 4656. The tail becomes 
very diffuse, and there is a color gradient along the ma- 
jor axis of the galaxy with the outermost regions of the 
tail becoming bluer. The exact shape and extent of the 
regions are shown in Figure [T] and detailed in Table [TJ 
while a color image of the system is shown in Figure [2j 
We note that in Figure [1] UV emission is visible further 
Southwest than the end of the tail region; this region is 
off the edge of the IRAC 3.6 and 5.8 /im images, so we do 
not include this area in all bands to preserve the colors of 
the spectral energy distributions (SEDs). For all bands 
and regions we adopt a detection threshold of S/N = 10. 

2.2.1. Surface gradients 

When computing photometry for sources near the de- 
tection limit, as necessitated here for most filters, small 
background fluctuations can add significant systematic 
errors to the results. These background gradients are in- 
dependent of the random errors used in our S/N thresh- 
olds and act to increase the true error in wavelength 
bands where an object has a small number of counts 
relative to the background level. 

To combat this problem we fit and subtract polyno- 
mial surfaces to all images before they are photome- 
tered. Generally when subtracting a background fit, we 
seek a solution that converges at large polynomial or- 
der. However, the commonly used scheme of masking 
sources tends to fail at high orders: restricting the lo- 
cation of points to perform the fit leads to poor sur- 
face subtraction in areas where a significant number of 
the pixels are masked, and increases errors overall when 
masks force non-optimal point selection (especially for 
polynomial surface fits). 

To solve these problems we took an iterative approach 
to build surface fits that were both accurate and stable 
at high polynomial orders. We started by fitting a zeroth 
order Lagrange polynomial (a constant offset) to the un- 
masked data using the IRAF task imsurfit. This fit is 
minimally sensitive to extended sources on the image. 
We replaced the pixel values of all masked regions by 
the value of the fitted zeroth order background and the 
resulting image is run through imsurfit with order = 1. 
This process is repeated while increasing the polynomial 
order until the change in the surface fit between orders 
asymptotes to a minimum value, which we find happens 
after 8 orders for most images (although for the 2MASS 
H-band we go to 19th order due to the larger variations 
in the sky background in that bandpass). 

For the GALEX, SDSS, and 2MASS images, the 
change between polynomial orders asymptotes at 
<0.05%. Because the IRAC and MIPS mosaics are not 
perfectly rectangular our algorithm's performance is not 
as good, with variations <1.5%, although we note that 



TABLE 1 

Apertures 



Name 


Region Type 


RA 


Declination 


Dimensions" 


Position Angle 






(J2000) 


(J2000) 


(arcseconds) 


(° 


E from N) 


NGC 4656UV 


Ellipse 


12:44:15.68 


+32:17:00.4 


180.3x115.8 




36.36 


NGC 4656 


Ellipse 


12:43:53.61 


+32:08:46.0 


390.0x142.2 




40.00 


NGC 4656-Hook 


Box 


12:44:08.82 


+32:12:36.8 


127.2x82.8 




100.0 


NGC 4656-Center 


Box 


12:43:57.30 


+32:09:47.6 


376.2x162.5 




30.00 


NGC 4656-Tail 


Box 


12:43:40.62 


+32:05:32.9 


328.5x202.8 




50.00 



For ellipses: the semi-major and semi-minor axes. For boxes: the length of the sides. 




Fig. 2.— False color RGB image of NGC 4656 and NGC 4656UV. Blue channel: GALEX FUV. Green channel: SDSS g. Red channel 
SDSS i. The SDSS images have been smoothed to a resolution of ^l". Red contours are for H I column density, drawn at 3, 8, 16, 32, 64 

cm '' 

its H I contours and those of its parent galaxy. 



128(7 where a = 2.61 X 10^^ cm ^. While both sides of NGC 4656 are extremely blue, only NGC 4656UV has distinct separation between 
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Fig. 3. — Demonstration of our surface-subtraction technique. 
Left column: original image. Right column: 9th order surface fit 
using our technique. From top to bottom the rows are: GALEX 
FUV, SDSS r, SDSS z, IRAC 3.6 ^m, and MIPS 70 fim. Green 
ellipses are the masking apertures used for the surface fits. All 
images have been rotated to North-up, East-left orientation and 
smoothed by a gaussian with (t=6 pixels to highlight faint struc- 
ture. Our algorithm accurately fits the surface gradients for all of 
the images at high polynomial order, even for images with signif- 
icant padding (e.g. the FUV, 3.6 ^m, and 70 fim images) and/or 
large mask regions (e.g. the r and z images). 

the largest variations are confined to the corners of these 
images and that the surface fit generally has less variation 
in the central parts of the mosaic. As Figure [3] shows, 
our method is robust regardless of the fitting order. 

While our method is a large improvement over non- 
iterative fits, there are still limitations to surface- 
removal. When the surface gradients are variable over 
scales comparable to the objects to be photometered, 
it is impossible to accurately subtract them without 
also subtracting the object. Fortunately in most of our 
images the gradients are on large scales and/or both 
NGC 4656UV and all parts of NGC 4656 have signifi- 
cantly higher counts than the background fluctuations. 
However, there are two exceptions to this which will be 
discussed further in i i3.4l and i i3.5l 

3. RESULTS 

We present our results from archival and survey obser- 
vations in two parts. First, we present the VLA 21 cm ob- 
servations and evidence from the cold neutral H I gas that 



NGC 4656/NGC 4656UV is a disturbed system. Then 
we present the UV/optical/infrared imaging and model- 
ing of the spectral energy distribution (SED) to demon- 
strate that the age and metallicity of the NGC 4656UV 
stellar system is consistent with the interaction history 
of the group, providing a self consistent picture for the 
creation of the tidal dwarf galaxy candidate. 

3.1. Evidence for the interaction history of NGC 4656 

The distribution and kinematics of H I gas in the en- 
vironment of the NGC 4656/NGC 4631 group provides 
the strongest evidence that the two edge-on spiral galax- 
ies have undergone a recent interaction. On the galaxy 
group scale, H I "fingers" extend from NGC 4631 towards 
NGC 4656, as if the H I gas has been stripped out of the 
former as the result of a gravitational encounter between 
the two galaxies. This has been well documented by two 
sets of authors, and we refe r the reader to Figures 1 and 
3 of l\yehachew et al.l (|1978[ ) and Figures 4 and 7 of lRandl 
()1994[) for a perspective of H I across the entire group. 

We use the archival VLA observations to investigate 
the kinematics and gas morphology of NGC 4656. We 
present the H I data cube in the form of velocity channel 
maps in Figures 2] and [SI total intensity map (moment 0) 
contours overlaid on a composite image in Figure [5J and 
an intensity weighted velocity (momen 1) map in Figure 
[6]). The moments are calculated independently for each 
pixel. The intensity weighted velocity map was created 
by clipping the emission below la of our RMS sensitivity 
in order to avoid large outliers. The systemic velocity of 
NGC 4656 is 645 ± 3 km s^\ while NGC 4656UV has a 
systemic velocity of about 570 kms~^. 

The channel maps show the presence of warping in 
the disk, as well as counterrotating and extraplanar gas. 
In each panel the "X" marks the center of the galaxy 
found through modeling of the rotation curve. In chan- 
nels 754.6 kms~^ through at least 723.6 kms~^ we see 
gas northeast of the center of the galaxy which is counter- 
rotating with respect to the disk. Channels 723.7 kms~^ 
through 703.0 kms^^ show warping in the outer portion 
of the disk. Channels 600.0 kms^^ through 548.5 kms^^ 
show evidence of extr aplanar gas (north of the "X" ; iden- 
tified as a "worm" bv lRandlll994[ ). The total H I intensity 
map (Figure ^ shows the neutral hydrogen gas extends 
beyond the optical plane of the galaxy both radially and 
vertically. In H I the TDG (red ellipse in Figures |3] and 
[S]) appears as a continuous extension of the galaxy to 
the northeast. We used our knowledge of the optical 
morphology to exclude the gas in this region from our 
modeling of the rotation curve of the parent galaxy. 

3.2. Modeling the rotation curve of NGC ^656 

We attempted to fit a rotation curve to NGC 4656 
using a t i lted ring model following the procedure of 
iBegemanI ()1989[ ) and iHess et al.l (|2009f ) with mixed suc- 
cess. The top panel of Figure [7] shows the results of 
fitting a rotation curve to concentric rings assuming an 
ideal, flat disk with constant inclination and position an- 
gle. We derived a systematic velocity of 645 ± 3 kms~^ 
and inclination of approximately 79° ± 3°. We see that 
the rotational velocity rises slowly with radius, reaching a 
maximum of ~ 65 kms^^ before roughly flattening out. 
From the rotation curve, we derived a total dynamical 
mass of 1.9 X 10^° M^ for NGC 4656. 
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Fig. 4. — Channel maps of NGC 4656 covering 754.6 km s~^ to 641.2 km s"'^ heliocentric radio velocity. Data are plotted at full spatial 
and spectral resolution. The beam size is drawn in the lower left corner of the first panel. Contours are plotted for 3, 10, 20, 40, and SOct 
where a = 5.03 X 10^® cm~^. Grayscale is in units of mjy beam~^. The "X" marks the center of the galaxy found through our modeling 
of the rotation curve. The red ellipse marks the position of the TDG candidate. Note the significant amount of gas to the northeast of the 
center in all panels, which is counterrotating with respect to the rest of the galaxy. 



The middle and bottom panels of Figure [7] demonstrate 
the quality of our assumption that the galaxy is an ideal 
disk by plotting how the inclination and position angle, 
respectively, vary with radius when we fix the rotation 
velocity of each ring to that found in the simple, ideal 
case. The ideal case fails to match all the features ob- 
served in the data: it cannot account for the twisted 
isovelocity contours in the intensity weighted velocity 
map (Figure ^ at the center of the main galaxy, indi- 
cating that NGC 4656 deviates from regular rotation. 
Further modeling shows the twisted isovelocity contours 
can be fit with a position angle which varies with radius, 
as demonstrated in the bottom panel of Figure El How- 
ever, when we examine the data cube, the center of the 
galaxy is dominated by non-circular motions and clumps 
of counterrotating H I gas at "forbidden" velocities, seen 
in velocity channel maps, which twist the isovelocity con- 
tours. 

Figure [5] shows a position- velocity (p-v) slice through 
the H I data cube along the major axis of the galaxy. 



This slice is atypical compared to what is usually ob- 
served in spiral galaxies, which e xhibit a qu ickly ris- 
ing, then flat rotation curve (e.g. iSofue fc R ubin 200l| 
and references therein). In the inner 300 arcseconds 
(~ 10 kpc) the gas exhibits a huge velocity dispersion 
(^ 100 — 150 kms~^) compared to the overall change in 
the rotation velocity with radius (of order 60 kms^^). 
The H I gas is not concentrated at the systemic velocity 
of the galaxy, but at higher and lower velocities as if the 
gas is expanding in the radial direction. Given the shal- 
low potential well inferred from the slowly rising rotation 
curve and that NGC 4656 is a fairly low mass spiral (100 
times less massive than the Milky Way), the modest star 
formation rate (a factor of 4-5 times less than the Milky 
Way, see iJ3.4.2p may be sufficient to blow gas out of the 
center of the galaxy. 

We believe the disturbed H I gas and gas at large 
heights above the disk of NGC 4656 are due to a com- 
bination of the galaxy-galaxy interaction and the cur- 
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Fig. 5. — As in Figure |4] but covering 630.9 km s~^ to 517.6 km s~^ heliocentric radio velocity. H I associated with the TDG is found 
most strongly in channels covering 600 km s~^ to 527.9 km s~^. Between this and Figure|4l it is clear the galaxy is rotating, but is very 
asymmetric. 



rent episode of vigorous star formation. Extraplanar and 
counterrotating gas has been observed in several edge-on 
and inclined galaxi es (e.g. NGC 891, [Q osterloo ct al. 
mrA NGC 240 3, iFraternali etHI 1200^ NGC 2997, 



Hess et al]l2009( ). The presence of this gas may be in- 



dicative of active star formation, however the largest 
H I clouds are likely evidence of accretion onto the 
galaxy — either of primo rdial origin (e.g. Hess et al. 2009; 
iWilcots fc Millenll99 8^. or as the result of tidal material 
falling back onto the dis k, such as has been modeled by 
IHibbard fc Mih^ (IT995I ). 

3.3. HI gas m NGC 4656UV 

Figure [2] shows an H I envelope spatially coincident 
with NGC 4656UV, although offset from the apparent 
UV-optical center of the TDG candidate. This H I en- 
velope is connected to its parent galaxy by a smooth ve- 
locity gradient in the H I gas: the p-v slice intersects the 
TDG candidate, identified as the density enhancement 
around —100 kms~^ (Figure [8]). 

We interpret the H I density enhancement and the ve- 



locity gradient as evidence of the H I gas having been 
pulled out of NGC 4656 during an interaction. The 
H I envelope serves as the gas reservoir out of which 
NGC 4656UV is currently forming stars. Summing the 
H I emission in this envelope yields a total fiux of 2.7±0.5 
Jy. At 7.2 Mpc this corresponds to a total H I mass of 
3.8 X lO'^ M0 for the TDG c andida te, consistent with the 
range of H I masses in the iDuc fc Mirabel (1999) sam- 
ple. The peak H I column density towards the TDG is 
1.3 X 10^^ cm~^, however, the column density coincident 
with star formation across the rest of TDG is below the 
supposed thr eshold of 10 M^d pc~^ for t he disk instabil- 
ity criterion ([Martin fc KennicuttI l200l|) . At the lowest 
level, we observe star formation occurring at gas surface 
densities of 1.6 Mq pc^^. It is important to note that 
Y^gas here is based only on H I observations, neglecting 
the molecular gas component which would raise the to- 
tal gas surface density. Nonetheless, the TDG candidate 
contributes to a growing body of evidence that star for- 
mation can occur in low density environments. 
The H I density enhancement we associate with 
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Fig. 6. — The intensity weighted velocity (moment 1) map of 
NGC 4656 and TDG candidate, NGC 4656UV. Emission has been 
clipped at 7cr, and contours range over 540 — 720 kms"^ in steps of 
10 kms"^. We note there is a continuous velocity gradient from the 
parent galaxy across the TDG companion. The twisted isovelocity 
contours show that the parent galaxy is clearly disturbed. The 
extraplanar emission to the northwest and southwest may be real 
emission associated wit h the H I "fingers" betw een NGC 4656 an d 
NGC 4631 reported bv lWeliachew et al.l HW!^ and lRandl JTOOi) . 

NG C 4656UV is evident in previous observations, how- 
ever iRandl ()1994[ ) proposed it was one side of an H I 
ring encircling the optically bright edge on galaxy. We 
consider this and other possible explanations in []4l 
The identification of the UV-bright /optically faint stellar 
component to the northeast of the galaxy (NGC 4656UV) 
is motivation for reconsidering the gaseous structure. 
There is no stellar counterpart to the southwest portion 
of the proposed H I ring. 

3.4. Spectral energy distribution of NGC 4^656UV: 
recent star formation and dust limits 

A summary of the UV, optical, and infrared photom- 
etry for NGC 4656UV, NGC 4656, and various morpho- 
logically distinct star forming regions within NGC 4656 
is shown in Table [H An upper limit indicates the flux 
that would be present from a source with a lOa detec- 
tion. The uncertainties are based on random errors only 
and do not include systematic effects (e.g. un-removed 
surface gradients) . Therefore our uncertainties should be 
treated as lower bounds on the true uncertainty. 

The z-band and 24 /xm images appear to be signifi- 
cantly affected by surface gradients. Photometry of the 
24 /xm image resulted in a better than lOcr detection 
for NGC 4656UV; however a visual inspection showed 
a structure more consistent with small-scale background 
variations seen elsewhere in the image than with the mor- 
phology of NGC 4656UV observed in other bands, so we 
do not use this "detection" in our analyses. There is 
clearly z band emission from NGC 4656UV in the SDSS 
image, but it appears to be super-imposed over a surface 
gradient with amplitude comparable to NGC 4656UV. 
We discuss this detection more thoroughly in > j3.4.1l For 
these two bands we also estimated the magnitude of the 
surface gradients on the image and report these as addi- 
tional errors on the integrated flux of NGC 4656UV in 
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Fig. 7. — Top panel: rotation curve for NGC 4656 derived from 
a tilted ring model. Middle panel: varying inclination angle with 
radius. Bottom panel: varying position angle with radius. The 
rotation curve is shallower than suggested by the velocity channel 
maps and does not capture the presence of anomalous or counter- 
rotating gas. The position angle declines by about 60 degrees from 
the center to the outskirts of the galaxy. This change in position 
angle is the model's attempt at fitting the twisted isovelocity con- 
tours apparent in the moment 1 intensity weighted velocity map. 
While the galaxy is disturbed at the center, the shape of these 
isovelocity contours is likely due to the countcrrotating gas instead 
of a warp in the disk. 
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3.4.1. Spectral energy distributions 



We plot the SED of NGC 4656UV with the SEDs of its 
parent galaxy in the left panel of Figure [HI To facilitate 
comparisons, we also show the SEDs scaled to have the 
same FUV flux in the right panel of Figure [HI All regions 
of NGC 4656 are actively forming stars, and have similar 
UV-optical colors to NGC 4656UV (with the exception of 
the z-band, discussed below). The SED of NGC 4656UV 
is very similar to that of NGC 4656 's tail, the other 
low-density star formation environment in this system. 
Both these systems have similar ratios of FUV-NUV-w 
band fluxes, and steep optical slopes. The hook also 
exhibits this morphology, although with a nearly non- 
existent Balmer break, in contrast to the strong breaks 
of NGC 4656UV and the tail of NGC 4656. The main 
region (and, due to the main region's signiflcantly higher 
brightness than the tail or the hook, NGC 4656 as a 
whole) has a shallower slope in both the UV and the op- 
tical, indicative of old stellar populations present in the 
center of the galaxy. All parts of NGC 4656 contain cold 
dust, although the tail appears to be deficient in 70 /Ltm 
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Fig. 8. — A position-velocity slice taken from the H I data cube along the major axis of the parent galaxy. The TDG candidate 
NGC 4656UV is the density enhancement to the left at approximately -0.65 kiloarcseconds. The parent galaxy appears to have some 
rotation, however the H I gas is clearly disturbed. We believe this is in part due to a previous gravitational interaction and to the vigorous 
star formation taking place in the parent galaxy. 



TABLE 2 
Integrated Photometry (mJy) 



Band 


NGC 4656 


Hook 


Center 


Tail 


NGC 4656UV 


FUV 


77.4 ±3.8 


24.9 ±1.2 


45.4 ±2.2 


9.89 ± 0.48 


3.14 ±0.15 


NUV 


77.1 ±2.1 


22.6 ±0.60 


48.8 ±1.3 


9.23 ±0.25 


3.22 ±0.086 


u 


91.0 ±0.12 


22.3 ±0.035 


63.2 ± 0.074 


10.1 ±0.064 


2.80 ±0.061 


g 


165 ± 0.072 


32.1 ±0.029 


122 ±0.057 


17.5 ± 0.030 


5.67 ±0.025 


r 


178 ± 0.094 


32.1 ±0.032 


137 ± 0.069 


17.3 ± 0.044 


5.56 ±0.041 


1 


172 ±0.14 


28.2 ±0.040 


139 ± 0.093 


14.9 ±0.067 


5.29 ±0.071 


z 


153 ± 0.46 


27.4 ±0.12 


134 ± 0.28 


6.81 ±0.28 


10.0 ± 0.26 ± 37'' 


J 


199 ± 8.2 


32.1 ±2.0 


159 ±4.9 


< 50.5 


< 58.8 


H 


206 ± 14 


< 34.6 


178 ± 8.3 


<85.6 


< 113 


K 


159 ± 13 


< 32.4 


126 ± 7.8 


<80.0 


<87.2 


3.6 


80.2 ± 0.027 


11.4 ±0.0079 


64.1 ±0.019 


6.38 ±0.014 


2.39 ±0.018 


4.5 


58.1 ±0.035 


8.65 ±0.0097 


46.4 ± 0.023 


4.10 ±0.019 


a 


5.8 


63.0 ±0.17 


9.22 ±0.050 


48.5 ±0.11 


10.0 ±0.10 


< 0.941 


8.0 


99.4 ±.17 


17.7 ±0.043 


70.1 ±0.10 


6.61 ±0.097 


a 


24 


530 ± .31 


135 ± 0.088 


387 ±0.19 


17.4 ±0.18 


9.51 ±0.17 ±23'' 


70 


9850 ± 100 


2290 ± 48 


6760 ± 82 


493 ± 23 


< 122 


160 


11800 ±170 


1890 ± 65 


841 ± 140 


669 ± 47 


<388 



'^ Significant portions of NGC 4656 were off the IRAC 4.5 and 8.0 /^m tiles. To avoid inac- 
curacies from making assumptions about the morphology of the un-observed regions (as 
well as issues with surface gradients along the edges of the images) we do not photometer 
NGC 4656UV in these bands. 

'' Significant small-scale surface gradients are co-spatial with the position of NGC 4656UV 
in the z and 24 fim images, which caused the random errors to be much smaller than the 
true errors. For these two detections we first report the random error, then an estimate 
of the systematic error based on the amplitude of other small-scale surface gradients in 
the images. 



emission compared to the rest of the galaxy. 

The z-band detection in NGC 4656UV is odd, con- 
sidering how similar the rest of its UV-optical SED 
is to the postulated parent galaxy. Significant excess 
flux in the z-band would be a strong indicator of an 
older stellar population, one either formed in-situ (if 
NGC 4656UV is a pre-existing dwarf galaxy) or ripped 
from NGC 4656 at the same time as the neutral gas cur- 
rently powering the ongoing episode of star formation 
(if NGC 4656UV is a TDG). A visual inspection of the 
z-band image clearly shows a signal that is consistent 
spatially with NGC 4656UV, however the detection is at 
small DN/pixel and there are similar-magnitude surface 
gradients nearby, similar to the 24 /im band (see row 3 
of Figure [31 and note especially the vertical bands on the 



image). We believe that these gradients are responsible 
for artificially boosting the z-band flux given the other- 
wise strong similarity between the SEDs of NGC 4656UV 
and NGC 4656, but in i i3.5l we consider the consequences 
of it being a real feature in NGC 4656UV's SED. 

3.4.2. Star formation rates 

We compute NGC 4656UV's star formation rate (SFR) 
using the conversion from FUV flux given by IKennicufrg 
(fT99l : 

SFR{Mq yr-^) = 1.4 x 10"^*i,.(erg s-^Hz-^). (1) 

At the assumed distance of 7.2 Mpc, NGC 4656UV's 
FUV luminosity is 1.95±0.09xl026 ergs'^ffz-^ There- 
fore, we calculate a star formation rate of 0.027 ± 0.001 
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Mq yr^^ for the TDG candidate. This is consistent 
with t he mean star formati on rate of TDGs accord- 
ing to iDuc fc Mirabell (|19990 . For the parent galaxy, 
NGC 4656, we compute a SFR of 0.666 ± 0.03 M© yr"!, 
roughly 25 times greater than that of NGC 4656UV. The 
FUV SFR for NGG 4656 is comparable to literature es- 
timates of the SFR from Ha, but roughly a factor of 
two smaller than the far- infrare d SFR and the 1.4 GHz 
radio SFR (jMapelh et al.ll2010t ). This discrepancy may 
be caused by dust attenuation on the UV e mission of 
NGC 4656, a source of uncertaint y no ted by iKennicutti 
(|1998[) : however, we will show in ^3.5\ that evolutionary 
synthesis models predict very low amounts of dust in 
NGC 4656UV, implying that attenuation corrections on 
the S FR of the T DG candidate should be small. 

The Kennicutt (1998) conversion formula has been op- 
timized for relatively stable, long-lived, quiescent galax- 
ies. Additionally, the calibration of Equation ([1]) depends 
on the initial mass function and stellar library chosen for 
modeling this relationship. As NGC 4656UV is a very 
different star-forming environment than the galaxies used 
to obtain Equation ([Ij the SFR computed using this for- 
mula should be treated with caution. 

3.4.3. Should we expect to see dust in NGC 4656UV? 

One of the most obvious differences between the SEDs 
of NGC 4656 and NGC 4656UV is the lack of dust de- 
tected from mid- and far-IR emission in NGC 4656UV. 
However, it is unclear whether the non-detections are due 
to a true lack of dust or the faintness of the TDG candi- 
date. Therefore, we attempt to estimate what we should 
expect to observe for NGC 4656UV based on its UV and 
optical properties. 

While all portions of NGC 4656 have similar SEDs 
in the UV and optical to NGC 4656UV, NGC 4656 's 
tail is the most alike. If we assume that this similarity 
extends into the mid-IR, then using the same ratio of 
NUV to 70 ^m flux as the tail of NGC 4656, we predict 
that NGC 4656UV should have a flux of 0.172 Jy at 70 
/xm: roughly 1.5 times the lOcr detection limit. We use 
the NUV flux in this ratio because both bands should 
trace the ongoing star formation event, however we note 
that no ratio of wavelengths predict a flux below our 
lOtT detection limit. The estimated NUV/ 160 /im flux 
ratios are, however, consistent with a non-detection for 
NGC 4656UV. These ratios indicate that NGC 4656UV 
has a lower cold dust density than even the most diffuse 
regions on NGC 4656. 

Given the very active star formation in NGC 4656UV, 
it is also valuable to estimate the expected flux from 
polycyclic aromatic hydrocarbons (PAHs). In a sam- 
ple of star forming r egions formed in collisional debris, 
iBoquien et al.l (|2009) report 8.0 /im detections in all sys- 
tems, indicating that many tidal features exhibit emis- 
sion from small grains of PAHs (although Holmberg IX, 
a nearby TDG candidate that will be discuss ed in more 
depth in i j4.2[ was not detected in 8.0 /im bv lDale et al.l 
120071) . PAH emission ha s also been detected in TDGs 
at several wavelengths bv lHigdon et al . (2006). Unfortu- 
nately, since the archival IRAC 8.0 fim image only con- 
tains about half of NGC 4656UV, a direct comparison 
of 8.0 /x m fluxes with UV em ission using the formulae 
given by IBoquien et al.l (|2009f) is impossible. However, 
the general lack of dust emission in the IRAC bands sug- 



gests an overall under abundance of PAHs. This may be 
due to NGC 4656UV's very low metallicity, consistent 
with our SED analysis and GALEV modeling (detailed 
in ^53)1 . 

3.5. Evolutionary synthesis modeling 

In order to gain more insight into NGC 4656UV's star- 
formation history, age, metallicity, and extinction we cre- 
ated a se t of stellar populati on synthesis models using 
GALEV (jKotuUa et all 120091 ). Being chiefly interested 
in the ongoing starburst we chose an exponentially de- 
clining burst star formation rate with an e-folding time of 
200 Myr. We considered a set of five metallicities rang- 
ing from [Fe/H] = —1.7 to 0.0 (solar), including a model 
with a "chemically consistent metallicity" . The "chemi- 
cally consistent" model starts with the burst at primor- 
dial metal abundance, then tracks the evolving gas-phase 
metal abundance as the stellar p opulations age , die, a nd 
return metals into the ISM (see iKotulla et al.l 120091 for 
more information). We include this model as the most 
limiting case of low metal enrichment. We varied the 
amount of dust, using a Calzetti attenuation model with 
extinction up to E(B-V) = 0.5 mag (Calzetti et al. 2000). 

For each model, fluxes are computed for all passbands, 
then a least-squares fit is performed at each model age 
to fit the model fluxes to the data for NGC 4656UV. 
To focus on the recent star formation we fit only to the 
\JY+ugri data, the region of NGC 4656UV's SED which 
would be dominated by a young, recently formed stel- 
lar population. In Figure [TU] we plot the model with 
the best fitting age for each metallicity, overlaid with the 
measured flux of NGC 4656UV in each wavelength band. 
Adding dust to the models decreased the resulting age of 
the system. For example, adding E(B-V) = 0.1 mag of 
dust resulted in a younger best fltting age for the TDG 
by a factor of 1.5 — 2, regardless of metallicity. However, 
the best fltting models for all metallicities had E(B-V) < 
0.05 mag, with most models requiring zero dust. As the 
metallicity increases there is a slight trend of increasing 
best flt age: the best flt ages are 280 Myr at [Fe/H] = 
-1.7, 292 Myr at [Fe/H] = -0.7, 320 Myr at [Fe/H] = 
—0.3, and 288 Myr at solar. The chemically consistent 
model has the youngest age, 264 Myr. The best fits over- 
all are given by the chemically consistent and [Fe/H] = 
— 1.7 models. At the best fit age the chemically consis- 
tent model has a metallicity of [Fe/H] = —2.36, ~0.65 
dex lower than even the least abundant fixed metallicity 
model. 

The apparent low metallicity of NGC 4656UV is sur- 
prising given that TDGs are traditionally metal-rich in 
relation to normal dwarfs. However, the system as a 
whole is metal-deficient with NGC 4656 having a metal- 
licity of -1.12 dex (Ma pelh et al.ll2010 ). NGC 4656UV's 
metallicity is ~10 times smaller than that of its par- 
ent galaxy. This may be explained if NGC 4656UV was 
formed out of gas stripped from the extreme outer disk 
of NGC 4656, which would be more metal-poor than gas 
closer to the center of the parent galaxy. While TDGs 
of all luminosities appear to have a f airly consistent 
Z0/3 metalhcity (Due & Mirabell 119990 . much higher 
than found for NGC 4656UV, the overall dearth of mat- 
ter and metals in this group indicate that NGC 4656UV 
is forming in a different environment than most TDGs. 

If we assume that the z-band excess is due to back- 



NGC 4656UV 



11 



3io-"» 



■"^t*^ 



-Tail 
■ Main 

oHook 



^* I ~H3+^^^ 



\ -\f*i. 



^ 



1 ' 


1 


1 






. NGC4656UV 




_ 




-Tail 


_ 






■ Main 


- 


- <-v 




oHook 


- 






^ ^^ 


















v^ / -. W 










M« 


"^ ^ ^ 




'V " -t^^ 


: ; >^ ^ ^ 






* \ ^ ^ 




' /' ^^ \ - 


\ ^ 


\ 


' ' ' ^ -. V - 


\\ ^ 


^ 


^11 *■ 


^\ 


\ \ 




\ 


x,^ uVl 


^^e^ , 




'er *^^^^ -' ' 1 


- 


^ 1 ^ / 


' 




\^^--- 


,-u 




Pen 




I 








1 10 

\{fj.m) 



X{fj,m) 



Fig. 9. — SEDs of all regions of NGC 4656 and NGC 4656UV. The left panel preserves the relative fluxes while the right panel scales 
all the SEDs to have identical FUV fluxes. Horizontal error bars denote the width of the filters. The entire NGC 4656-4656UV system is 
clearly undergoing a siinilar-agc burst of star-formation. 
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Fig. 10.— GALEV fits to the bluest six bands of NGC 4656UV's SED. Data from NGC 4656UV are shown as black points, while the 
best fitting GALEV models to the UV+ugri data for each metallicity bin are shown as colored lines. Blue, green, and purple lines denote 
the best fitting models, while the orange line is for the most metal-poor model at an age of 13.6 Gyr. The most metal poor model flts to the 
UV and optical data well reproduce both the UV and optical data but are also consistent with the 3.6 /im detection. However, attempting 
to flt the z-band flux by adding a highly evolved stellar population fails, ovcrpredicting the flux in neighboring bandpasses. The red line 
denotes the spectrum of an M giant template star from iBruzual fc Chariots 1^199^), which would be the characteristic spectrum of any old 
stellar population that could fit the z-band fiux without overpredicting r and i band emission. 
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ground gradients we find that while the models are min- 
imized only to the bluest six detections, both the chem- 
ically consistent and [Fe/H] = —1.7 models do a fairly 
good job approximating the 3.6 ^m flux. The chemically 
consistent model indicate that all of the 3.6 /im emis- 
sion comes from starlight while the [Fe/H] = —1.7 model 
allows for a dust contribution of roughly 10-25% of the 
total light. At [Fe/H] = -0.7 and higher the models 
over predict the 3.6 //m emission, a further indication 
that NGC 4656UV contains low metallicity gas. 

We investigate the possibility that the z-band excess 
is due to an older, more evolved population from a much 
earlier starburst by over plotting the SED of the [Fe/H] — 
-1.7 model for a 13.6 Gyr old burst on FigurefTOl Even for 
this extremely old population, too much light is emitted 
in the r, i, and 3.6 /im bands to fit NGC 4656UV's SED. 
In order to fit the z-band signal without over predict- 
ing the emission in other bandpasses the old population 
would have to have the characteristic spectrum of a late 
M-type star (plotted in red in Figure [TOl) . The inability 
of the 13.6 Gyr model to even approach a good fit of the 
data is another indicator that the z-band flux is aber- 
rant. However, we note that even with the current data, 
we cannot conclusively rule out the possibility that there 
is an extremely old stellar population in NGC 4656UV. 

4. DISCUSSION 

While it is generally difficult to confirm that a galaxy is 
indeed a tidal dwarf witho ut detailed stellar population 
studies (e.g. Holmberg IX, iSabbi et al."2008') or detailed 
dynamical modeling (e.g. Hibbard & Mihos 1995), we 
favor NGC 4656UV's recent formation through a tidal in- 
teraction, most likely between NGC 4656 and NGC 4631. 
In the following sections we discuss the likelihood of other 
possible origin scenarios in the context of the kinematic 
and stellar data we have presented. Finally, we compare 
NGC 4656UV to observations of other star forming sys- 
tems and tidal dwarf candidates. 

4.1. NGC 4656UV: the case for a tidal dwarf 

We postulate that NGC 4656UV formed out of gas 
stripped from the outskirts of NGC 4656 during a recent 
interaction. The most likely culprit for this interaction is 
NGC 4631. The difference in the line-of-sight velocities 
towards NGC 4656 and NGC 4631 is -40 kms^^, which 
is extremely low for a group velocity dispersion, and may 
suggest that most of the relative velocity between the 
two galaxies is in the plane of the sky. With a projected 
separation of 36 kpc and assuming a typical group veloc- 
ity dispersion of 200 kms""'^, we estimate that the two 
spiral galaxies could have interacted as recently as 230 
Myr ago, comparable to the post-starburst age estimates 
from the GALEV models. Detailed kinematic modeling 
is required to strengthen the case that NGC 4656UV was 
created by an interactio n, b ut the H I br idge observed by 
iWeliachew et all (|1978f ) and|Rand| (|1994r ) and lend weight 
to NGC 4631 as the perturber. 

Several features of the UV-IR SED and the neutral 
hydrogen maps imply a tidal origin for NGC 4656UV. 
The disassociated morphology of NGC 4656UV, seen 
both in H I contours and especially in the UV stellar 
emission, indicates that NGC 4656UV is a separate en- 
tity from NGC 4656. NGC 4656's tail, while possessing 



a similar SED to NGC 4656UV, shows no signs of be- 
ing disconnected from the main body of the galaxy. At 
the same time, the existence of a light bridge between 
NGC 4656UV and its putative parent as well as the con- 
tinuity of the rotation of the neutral hydrogen gas be- 
tween the systems implies a strong dynamical connection 
between them. 

If we assume the velocity gradient across the H I den- 
sity enhancement is due to rotation in the TDG candi- 
date, we can use the measured baryonic mass, velocity, 
and size to estimate the virial ratio 



T _ RV^ 



(2) 



where T and U are the kinetic and potential energies, 
respectively, and we have made the simplifying assump- 
tion that NGC 4656UV is a constant-density sphere. We 
compute a total H I mass of 3.8 x 10^ Mq for the TDG. 
Using the g-band magnitude, the estimated distance to 
the TDG, and a M/L ratio of 1, we estimate a stellar 
mass of — 9 X 10^ Mq, giving a total H I -I- stellar mass 
of 4.7 X 10^ Mq. Assuming a diameter of 296" (~ 10.3 
kpc), equal to the average of the semi-major and semi- 
minor axes of our photometric aperture (large enough 
to contain most of the H I emission), and a total veloc- 
ity gradient of 73 kms~^, we estimate a virial ratio for 
NGC 4656UV of 1.13. This indicates that NGC 4656UV 
is right at the cutoff of being gravitationally unbound, 
perhaps unsurprising given that it appears to have been 
so recently liberated from NGC 4656 but has remained 
as a single entity. Combined with our estimate of the 
virial ratio, the large diameter of NGC 4656UV may in- 
dicate that the system is at least partially (if not fully) 
unbound. We note that this estimate is extremely rough: 
the density is not constant across the entire TDG, and 
the peak of the density distribution is not coincident 
with the center of the stellar distribution. This estimate 
also does not take into account kinetic energy due to a 
velocity dispersion, or expansion/contraction of the gas 
(jBournaud et al.ll2004l ). Additionally, even if some or all 
of NGC 4656UV is self-gravitating it is unclear whether 
or not this system will be fully ejected from NGC 4656 
or will fall back onto its parent galaxy. 

We can also use these estimates of rotation velocity 
and size to estimate a total dynamical mass: 



M, 



dyn 



G 



(3) 



We estimate that IvUyn « 1.6 x 10^ Mq for NGC 4656UV, 
only 3.4 times larger than the mass of the H I gas -I- 
stars. This discrepancy is similar to that found in other 
syste ms believed to be formed out of coUisi onal debris 
(e.g. iDuc et all [20071 : iBournaud et al.ll2007| ). Some of 
the missing mass may be accounted for by the presence 
of molecular gas (e.g. CO, H2). Regardless, it appears 
that NGC 4656UV has a relatively low ratio of dark-to- 
luminous matter compared to other dwarf galaxies (see 
Ashman 19_92_ for a review) , further supporting our hy- 
pothesis that NGC 4656UV formed from tidal interac- 
tions. 

Unfortunately, we cannot make a confident statement 
about the future of this star forming system and whether 
it might be long lived. Although the data suggests that 



NGC 4656UV 



13 



■V 

3 

DO 

a 



oTidal Tails 
■Bridges 
aTDGs 
•Other 



i 



Arp 
Arp 242 



I i-Sh"? s ■ 
m 



ii4fe'- 



%h^^ 



• '*' ^ArplSI 



NGC 4656UVA 



{Arp 2 



-0.5 0.5 1 

g-r (AB magnitudes) 

Fig. 11. — FUV-g v s. g-r color-color plot for all tidal features 
in ISmith et"all l(2010h . NGC 4656UV and Holmberg IX are also 
shown. Features identified by Smith ct al. (2010) as tidal tails 
are yellow open boxes, bridges arc red filled boxes, TDG candi- 
dates are displayed as blue open triangles (with the exception of 
NGC 4656UV, which is a filled green triangle), and all other objects 
are black circles. 

the system is not dark matter dominated, CO data would 
be even more powerful in revealing the internal kinemat- 
ics of the syste m and determining whether it is gravita- 
tionally bound (jBraine et al.ll200ll ). 

4.2. Comparison to other TDGs 

It is clear from the photometry and H I data that 
NGC 4656UV is vigorously turning gas into stars. 
We quantify this extreme star formation by comparing 
NGC 4656UV with a relevant sample from the literature. 
ISmith et all ()201CI( ) selected 42 nearby pre-merger, in- 
teracting galaxy pairs and obtained broadband GALEX 
photometry for the galaxies as well as their tidal tails, 
bridges, and likely tidal dwarf candidates. Many of these 
features are very bright in the UV and nearly invisible 
in the optical (e.g. see Arp 72 and 82 in Fi gure 2, Arp 
202 i n Figure 9, and Arp 305 in Figure 14 of .Smith et al.l 
l2OT0h . 

In Figure [Til we p resent the g-r vs. FUV-g colors of 



all lSmith et al.l ()2010f ) features with SDSS data compared 
to NGC 4656UV. We include Holmberg IX, a similarly 
young TDG near M81, for whose designation is stronger 
than most candi date TDGs due to resolved stellar pop- 
ulation studies ()Sabbi et al.l 120081) . Since no published 
SDSS photometry of Holmberg IX exists in the litera- 
ture, we photometered Holmberg IX in SDSS g- and r- 
bands, and obtained AB magnitudes of 15.18 ± 0.08 in g 
and 15.25±0.05 in r. We applied Galactic extinction cor- 
rections based on values given in NED*, and error bars 
in Figure ITT] include the uncertainty in transforming the 
corrections into g and r magnitudes. 

Even among tidal features that exhibit properties of 
very young stellar populations, NGC 4656UV is one of 

4 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 



the bluest objects. A visual inspection shows that the 
other very blue tidal dwarf candidates — the TDGs in 
Arp 202 and Arp 305, as well as Holmberg IX — are 
similar in appearance to NGC 4656UV. They all ap- 
pear clearly in GALEX, but are nearly invisible in SDSS. 
The apparent ^1 magnitude gap between NGC 4656UV, 
Holmberg IX, Arp 305, Arp 202, and the rest of the 
ISmith et al.l (|2010l) sample is intriguing, although with- 
out a larger sample it is impossible to tell if it is signifi- 
cant. (Note that NGC 4656 's tail would bridge this gap 
if plotted on Figure ITTj) . The gap between the bluest and 
reddest TDG candidates may indicate that the red TDGs 
are a more evolved form of NGC 4656UV-like objects, or 
that there are two separate populations of TDGs. 

As the only other known candidate TDG within 10 
Mpc, comparisons between NGC 4656UV and Holm- 
berg IX are especially interesting. From resolved stel- 
lar population studies, Holmberg IX boasts a very young 
stellar population, directly c omparable to the estimate 
of - 270 Myr for NGC 4656 (ISabbi et al.l l 2008l). and un- 
like m any of the "classical" TDG candidates ()Duc et al.l 
I2000f ). neither are located at the end of an o bvious, long 
(sever al 10s of kpc) tidal tail. Although ISabbi et al.l 



2008) could not completely rule out the presence of 
a skeleton red giant branch population of stars, Holm- 
berg IX is one of the most well-established TDG can- 
didates. The 10 Mpc threshold is roughly the limiting 
distance for which HST can do resolved stellar popula- 
tion studies, making NGC 4656UV one of the few known 
tidal dwarf galaxy systems eligible for follow-up HST ob- 
servations. 

4.3. Other possible origin scenarios 

We cannot rule out the possibility that NGC 4656UV 
is a pre-existing, low surface-brightness, gas-rich dwarf 
galaxy that has been highly rejuvenated by its inter- 
action with NGC 4656, or by the interaction between 
group members NGC 4656 and NGC 4631. This scenario 
can be tested attempting to detect an old stellar popula- 
tion's red giant stars through res olved stellar popula tion 
studies (similar to that done by ISabbi et al.l ()20Q8f ) on 
Holmberg IX), by determining NG C 4656UV's molecu- 
lar gas content (JBraine et al.ll200ll ). or by conclusively 
determining its metallicity throu gh, for example , oxy- 
gen abundances in H II regions (jDuc fc Mirabell 11994 
Weilbache r et al]l2000l ) . A confirmation of the low metal- 
licity estimated from the evolutionary synthesis model- 
ing could be an argument against NGC 4656UV forming 
strictly from tidal debris of its parent galaxy. 

We consider that NGC 4656UV is actually a pro- 
jected tidal tail, viewed at an angle such that it ap- 
pears to be a real overdensity of stars and gas. We 
count this scenario as unlikely based on the velocity gra- 
dient of NGC 4656 /N GC 4656UV shown in Figure El 
iBournaud et al.l (|2004D use simulations to show the sig- 
nature of long tidal tails viewed edge-on in projection is 
revealed by a change in the velocity gradient along the 
line of sight: the velocity appears to rise with radius to 
a point before falling back towards the systemic veloc- 
ity of the parent galaxy. Although we cannot completely 
rule out this scenario, the p-v diagram (Figure [S]) shows 
that the velocity across the major axis of the TDG can- 
didate may be continuing to increase with distance from 
the center of NGC 4656. 
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Finally, we cannot rule out that NGC 4656 is a UV- 
bright stellar complex in the outer disk of NGC 4656. 
These regions, uncovered by GALEX observations, are 
proving to be common in low density en vironments in 
the extreme outer disks of spiral galaxies (iThilker et alj 
[2005, 2007]) . They reside beyond the classic optical star- 
forming disk, and trace the extended spiral arm struc- 
ture revealed by H I observations. These stellar com- 
plexes are found in the peaks of the H I distribution, 
however, the total gas surface density is still below the 
classic star fo rmation threshold traced by the Ha and 
Tfgas relation (jMartin fc Kennicuttl 120011 ). In fact, these 
regions of star formation frequently do not have asso- 
ciated Ha emission. Similarly, there is no obvious Ha 
emission a ssociated with NGC 4656UV (see Figure 1 of 
iDonahuee t al. 1995). However, the absence of Ha could 
be due to the fact that young stars from a burst of star 
formation are visible for a longer period in the UV. 

4.4. Importance of UV data for TDG detection 

UV-bright TDGs like NGC 4656UV demonstrate that, 
in order to make a complete census of both tidal de- 
bris around interacting systems and the star formation 
density in the local Universe, the ultraviolet part of the 
spectrum cannot be ignored. A UV-selected sample is 
needed to further our understanding of these enigmatic 
dwarfs which are too faint in the SDSS wavebands to be 
easily detected in an optical survey. Four of the 10 TDG 
candidates in Figure [11] with UV and optical photom- 
etry lie in a different area of color-color space than all 
other tidal features, with FUV— 5 w 5 — r » 0. While it 
is unclear exactly what fraction of all TDG candidates 
are this blue, the relatively large number of UV-bright 
TDG candidates in this small sample imply that this 
population is not neglig ible. TDGs may be common: 
iHunsberger et al.l ()1996l ) estimate that up to half of all 
compa ct group dwarfs are TDGs, while [Due fc Mirabell 
()1999( ) indicate that a significant fraction of dwarfs in 
clusters could be recycled objects. Having an accurate 
census of the byproducts of galaxy interactions is criti- 
cal to understanding the evolution of groups and clusters, 
but without an unbiased UV survey a significant fraction 
of these tidal remnants will remain hidden from view. 

5. CONCLUSION 

We have identified a candidate tidal dwarf galaxy 
around NGC 4656 in GALEX near- and far-UV images. 
A significant amount of multi- wavelength archival data 
covering this object exists, including images from the 
FUV to 160 fixn and an H I data cube from the VLA. 
NGC 4656UV is chiefly characterized by an asymmetric, 
spatially distinct H I envelope, an extremely blue SED, 
and a lack of infrared emission associated with dust. 
The TDG candidate is currently undergoing a significant 
burst of star formation, making it extremely bright in the 
UV. We are unable to detect any emission in survey and 
archival data further red than 3.6 ^m, suggesting that 
there is little to no detectable large grain or PAH dust 
emission in this system, although deeper observations are 
warranted. 

The H I morphology of the entire NGC 4656 system 
is extremely disturbed featuring warps as well as coun- 
terrotating and extraplanar gas. We attempted to fit a 
rotation curve to NGC 4656 and found that this galaxy 



does not have the regular rotation expected from nor- 
mal spiral galaxies. The HI gas is extends high above 
the disk of NGC 4656 and has a large velocity disper- 
sion, suggesting a previous gravitational encounter and 
a strong recent burst of star formation. 

We used the H I velocity gradient across NGC 4656UV 
to make an estimate of the dynamical mass and the virial 
ratio, finding that for a set of simple assumptions the 
TDG may be either bound or unbound and that it has a 
low ratio of dark-to-luminous matter. Both of these re- 
sults indicate a system that has been formed via galaxy 
interactions and not from its own local matter over 
density. By fitting evolutionary synthesis models from 
GALEV to the spectral energy distribution of the tidal 
dwarf galaxy candidate, we find that NGC 4656UV's 
ongoing starburst is best represented by an extremely 
low-metallicity model with an age of '^270 Myr and no 
dust extinction. Supported by a simple dynamical argu- 
ment, we propose that the gas out of which NGC 4656UV 
is forming stars was stripped from its parent galaxy, 
NGC 4656, in a tidal interaction with fellow group mem- 
ber NGC 4631 approximately 200-300 Myr ago. 

When compared against other tidal features of inter- 
acting galaxies, we find that NGC 4656UV is bluer than 
tidal ta ils and bridges in a relatively large sample of tidal 
debris ([Smith et al.ll2010l) and it lies in a distinct region 
of the plot with three other TDG candidates. It is un- 
clear whether these blue candidate TDGs represent the 
youngest tidal dwarfs or are a distinct population from 
TDG candidates with redder colors. 

NGC 4656UV is an intriguing example of star for- 
mation in a low density environment, revealed through 
observations spanning almost the entire electromagnetic 
spectrum. Regardless of whether they are primorial in 
origin or formed from second-hand material, diffuse, gas 
rich systems like NGC 4656UV may be important sites 
of ongoing star formation in the nearby Universe. 
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